Mesoscopic Coulomb drag, broken detailed balance and fluctuation relations 
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When a biased conductor is put in proximity with an unbiased conductor a drag current can be 
induced in the absence of detailed balance. This is known as the Coulomb drag effect. However, even 
in this situation far away from equilibrium where detailed balance is explicitly broken, theory predicts 
that fluctuation relations are satisfied. This surprising effect has, to date, not been confirmed 
experimentally. Here we propose a system consisting of a capacitively coupled double quantum dot 
where the nonlinear fluctuation relations are verified in the absence of detailed balance. 
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Introduction — Mesoscopic physics offers a unique lab- 
oratory to investigate the extension of equilibrium- 
fluctuation dissipation theorems into the non-linear non- 
equilibrium regime The equilibrium fluctuation- 
dissipation theorem and its closely related Onsager sym- 
metry relations 0] are a corner stone of linear transport. 
It has therefore been natural to ask whether such rela- 
tions exist also if the system is driven out of the linear 
transport regime. For steady state transport, fluctuation 
relations have been developed which relate higher order 
response functions to fluctuation properties of the sys- 
tem P, • For example the current response to second 
order in the voltage (the second order conductance) is 
related to the voltage derivative of the noise of the sys- 
tem and, in the presence of a magnetic field, to the third 
cumulant of the current fluctuations at equilibrium [l], [j| . 

Clearly tests of non-equilibrium fluctuation relations 
are of fundamental interest. From a theoretical point of 
view, the task is to propose tests in which crucial rela- 
tions valid at equilibrium fail in the non-linear regime 
and to demonstrate that, despite such a failure, fluctua- 
tion relations hold. For instance we have suggested ex- 
periments which test fluctuation relations for systems in 
the presence of a magnetic field and in a regime where the 
Onsager relations are already known to fail Q], 0] ■ Such 
experimental tests are currently under way [8(. Here we 
propose to test fluctuation relations in a system where 
away from equilibrium we have no detailed balance. We 
consider two quantum dots in close proximity to each 
other such that they interact via long range Coulomb 
forces. The absence of detailed balance is manifest in a 
Coulomb drag [9j : the charge noise of one of the systems 
(the driver) drives a current through the other unbiased 
system [k|. Therefore, the drag current is a direct indi- 
cation that this fundamental symmetry is absent. Nev- 
ertheless, we below demonstrate that there exist fluctu- 
ation relations. 

The interaction of two systems in close proximity to 
each other plays a role in many important set-ups in 
physics. We recall here only the interaction of a detector 
with a system to be measured 11| which also provides a 




FIG. 1: (Color online) Sketch of two capacitively coupled 
quantum dots each one attached to two different terminals. 
For very large intradot charging energy, only four charge 
states are allowed, as depicted. Their dynamics is governed 
by the tunneling rates Tf 1 and 7; . 



test of fluctuation relations 12]. The shot noise current 



current correlations in nearby quantum dots has been 
measured by McClure et al. [13[ and discussed theoreti- 
cally [13, EH . Recently, reciprocity relations of two cou- 
pled conductors were proposed by Astumian jTij]. Here 
we emphasize that one conductor even if unbiased can 
act as a gate to the other conductor. As a consequence, 
the currents are not a function only of voltage differ- 
ences applied to each conductor but also depend on po- 
tential differences of one conductor to the other one. In 
an instructive work, Levchenko and Kamenev discuss the 
mesoscopic Coulomb drag for two quantum point con- 
tacts in close proximity [17] . In this geometry, charging 
of the point contacts can be neglected and the coupling 
of the two conductors is extrinsic via the capacitance of 
the leads. 

General theory. — The probability P(N,t) that N = 
(Ni, . . . , Nm) particles are transmitted through M leads 
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during time t charaterizes the statistical properties of 
our system. It is useful to consider the generating func- 
tion which is the logarithm of the "Fourier transform" 
T( X ) = IkJ2n P{N,t)e^ N of the distribution function. 
Here X is the vector of the counting fields. From the 
generating function all cumulants can be obtained by 
simple differentiation. The fluctuation relations are a 
consequence of symmetries of the generating function 
[H, IH ■ In particular (in the absence of a magentic field) 
it holds F{ix) = F{—iX + oV/kT), which is equivalent 
to P(N) = e NqV ' kT P{-N). Here qV/kT is the affin- 
ity vector with components given by the applied voltages 
Vi. By expanding the current through lead i, Li = (Jj), 
where li is the current operator, and the current corre- 
lations Sij = (AijATj), where Ali = ij — ij, 
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we can relate linear response current and equilibrium 
fluctuations by means of the fluctuation-dissipation the- 
orem, S 1 ?? = 2kTGi } j. The generalization to the weakly 
nonlinear regime reads 0, US EH 

S a p,*t + S a -y,p + Sfij.a — kT(G a ^i3j + Gp^ ai + G 7iQ ^).(3) 

Notably, we find that these nonlinear fluctuation rela- 
tions are valid even in the absence of detailed balance. 

To determine the general current-voltage character- 
istics and the nonlinear fluctuations relations for two 
interacting conductors, we employ the classical treat- 
ment of the Coulomb interaction that respects charge 
conservation (gauge invariance). We take the inter- 
action to be intrinsic, determined by the charges on 
the mesoscopic conductors, and assume the leads to 
be metals with perfect screening. Then, the dynamics 
of the system is determined by the sequential tunnel- 
ing between states with a well defined charge occupa- 
tion which obey the master equation p{t) = A4p(t) for 
the occupation probabilities 19]. Analogously, one can 
write the equation of motion for the generating function, 
C(<) - E N P( N ^y Nx ^ g iv en by ((t) = M( X )((t). The 
cumulant generating function T is given by the eigen- 
value of Jyt(x) that develops adiabatically from zero with 
small x 120] • Generally, the explicit expression for T 
is difficult to handle, so in practice it is more conve- 
nient to calculate the cumulants recursively order by or- 
der [llKi^. Thus, from coefficients cn k \ of the expansion 

T = J2{i k } c {h} ( e ' Xl ~ I)' 1 • • • (e lXM -l) lM , we obtain the 
current-current correlations up to any order 2l| : e.g., 
the current, I t = qJ2{i k } c {i k } S h^T.i k s, the zero fre- 
quency noise, Su = qh + 2q 2 J2{i k } c {l h } S h,2^'£,h,2, and 
the cross-correlations Sij = q 2 J2{i k } c {h} S h^i j ^T.h,2, 
where {Ik} = {h ■ ■ ■ Im} and q is the electron charge. 



Drag current and fluctuation relations. — In the follow- 
ing we explicitly show, using the previous formalism, 
the fulfillment of the fluctuation relations in a nonequi- 
librium system where detailed balance is broken. We 
consider two capacitively coupled two-terminal quantum 
dots (see Fig. Q]) with large intradot charging energy. 
The interdot coupling is described with a capacitance 
C . Hence, the dynamics is characterized by four charge 
states: the empty state |0) = 1 00) , the singly occupied 
states |u) = 1 10} and |d) = 1 01 ) , and the doubly occupied 
state |2) = 1 1 1 ) . Quite generally, the tunneling ampli- 
tudes are energy dependent. Therefore, we distinguish 
Ti, which denotes a tunneling process through barrier 
I = 1 , . . . , 4 when the system is empty, and 7; , which cor- 
responds to tunneling when the coupled dot is already oc- 
cupied. This is the minimal charge model that manifests 
violation of detailed balance leading to drag currents. 
Detailed balance is broken when the probability to trans- 
fer one charge from left to right differs from the reverse 
process (from right to left). For instance, an electron 
is transported from left to right in the drag system by 
the sequence |0)— »|u)— >\2)— >\d)— >\0) with a probability 
oc Ti72 whereas the probability to transport it from right 
to left is cx 71 F2. Clearly, both probabilities differ and a 
nontrivial current, the drag current /drag oc — 71 T2, 
will be generated. We need that (i) both empty and 
doubly occupied states are taken into account and (ii) 
the tunneling rates depend on the charge state. Thus, 
a model with three charge states only (|u), |d) and |0) 
or 1 2)) cannot break the detailed balance and the drag 
effect is absent. The biased dot then acts merely as a 
fluctuating gate on the other dot. 

For the system depicted in Fig [TJ writing £ = 
(Co, Cuj Cd, C2), the equation ( = M(x)C becomes, 
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where f± = E^e^lf, and 7 ± = E« ea e ±<x "yf, 
u = {1,2} and d = {3,4}. The tunneling rates m(-) 
and out(+) of the dot read Tf 1 — Tif^ and — JifjJ 
with /± = /±(w„ - qVi) (n = 0,1). Here, f+(e) = 
1 — /(e) and / _ (e) = /(e) denote the hole and elec- 
tron Fermi functions, respectively. The effective level of 
dot a with bare level e Q when dot /3 ^ a is uncharged 
(n = 0) is n a0 = e a + [q 2 C^ a /2 + q(Cxp Eie« C l V l + 
CJ2iep CiVi)]/CC, where Ci is the capacitance of the Zth 
barrier, C Eq = J2 lea d + C and C = (C Su C Sd - C 2 )/C. 
In the charged case (n — 1), we find fJ, a i = /i Q o + ^C with 
Ec = 2q 2 /C the energy needed to add a second electron. 

We now investigate the drag current, for which we take 
the upper subsystem as the drag circuit (Vi = V2) and 



the lower one as the driver. Then, I\ 



-h = ^dra R and 
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FIG. 2: (Color online) Voltage dependence of (a) the drag 
current /drag through the upper dot (Vi = V2) (b) the drive 
current /drive for the lower dot, (c) the current-current cor- 
relations in the drive system, S22 = (A/2A/2) and (d) the 
cross-correlation between currents at different conductors, 
524 = (A/2A/4), for the drag configuration, Vi = V2. Pa- 
rameters: Ti — = r, except 71 = O.ir, kT = 5hT, 
q 2 /d = 20W\ q 2 /C = 50W\ e u = e d = 0. 



we find 
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± (5) 
Ai ± 



(r 3 /£ +r 4 /±)/r d and = (73/3=1 +7*/«)/7d are 
nonequilibrium distribution functions. 

When the drive voltage V3 — V4 is small, detailed bal- 
ance must be broken also in the drive circuit in order to 
have a linear / d ra g : G 2 a °c (7i r 2 - 72ri)(74r 3 - 73P4). 
Therefore, asymmetry in both the drag and the drive sys- 
tems is required for a nonzero linear drag current. More- 



G^l, satisfying the Onsager-Casimir 



over, we get G 2 X \ 
reciprocity relations [2j. Note that if the drive conduc- 
tor is also unbiased (V3 = V4), equilibrium fluctuations 
are expectedly not enough to induce a net current. This 
can be seen in Fig. EJa). For low voltages there is a 
Coulomb gap where transport is not allowed. This result 
also demonstrates that the voltage difference between the 
two subsystems V\ — V3 plays a crucial role, affecting the 
dynamics: In this case, one of the conductors acts as a 
gate on the other one. The gate effect of the drag circuit 
onto the driver is shown in Fig. EJb), where we obtain a 
typical Coulomb blockade stability diagram for the drive 
current J 3 = — 7 4 = 7 drive . 

It is worth noticing that, at high enough drive bias, 



/drag is suppressed since the interdot capacitance brings 
the dot states outside the transport window. Then, the 
drag current peaks at an optimal value of V\ — V3 and van- 
ishes away from it. On the other hand, at very low tem- 
perature /drag is finite only within a voltage range defined 
by a*io < qVi < /in and min{gF 3 , qV A } < ^10,^11 < 
max{qV3, qVi}. As expected, the drag current increases 
with C, but the voltage window where /drag is observable 
becomes narrower. Then, for large coupling the drag cir- 



cuit effectively induces dynamical channel blockade [23 1 
in the driver and, eventually, the drive current shows 
electron bunching. 

If the drive system is symmetric, the sign of /drag de- 
pends on the asymmetry factor (I172 — 71^) due to the 
competition of processes transferring an electron in each 
direction, independently of the direction of /drive- These 
two contributions have been detected separately in cou- 
pled double dot systems in the cotunneling regime giving 
rise to bidirectional drag j24[. Note that the asymmetry 
of the drag system can be enough to get a negative drag. 

We now investigate the nonlinear fluctuation relations 
for our system. We first analyze the occurrence of /drag 
and the current cross-correlations Sij for different con- 
ductors (e.g., i — {1,2} and j = {3,4}). The observa- 
tion of drag current in one conductor requires the oc- 
currence of correlated tunneling events between the two 
dots involving the states |0) and |2). These correlated 
events lead to finite cross-correlations. This would not 
be the case for a model that includes only three charge 
states. Our minimal model of four charge states does 
generate correlations between the currents through the 
two dots. For example, at equilibrium, the fluctuation- 
dissipation theorem relates the linear drag current to 
the equilibrium cross-correlations for different conduc- 
tors, G2,4 = S^/2kT. Similarly to /drag, if both conduc- 
tors are symmetric, i.e., 71 T2 — 72 Ii and 73 E4 = 74 P3, 
Sij vanishes to first order in a voltage expansion. Figure 
H£d) shows that the cross-correlation between the drag 
and drive currents is finite only when there is a drag cur- 
rent flowing in the upper conductor. In general, the sign 
of the cross-correlations is not determined by the direc- 
tion of the averaged currents However, in our case, 
the cross-correlations are positive whenever the two cur- 
rents flow in the same direction, and negative when they 
are opposite. Interestingly, /drag can present negative ex- 
cess noise, i.e., the noise S22 decreases in the presence of 
drag, as shown in Fig. [2{c). S22 reaches its maximal 
value when the effective upper dot level is aligned with 
the Fermi level (25|. 

Finally, we explicitly check that these fluctuation re- 
lations [l|, 0-11] hold even for our system in which de- 
tailed balance is violated. Charge conservation in each 
subsystem implies I a = — / a and S aa = Saa = —S a a 
for two different terminals in the same conductor. Then, 
from Eq. ([3]) we derive the nonlinear fluctuation relations 
involving terminals of the same conductor, and rewrite 
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FIG. 3: (Color online) Verification of the non-equlibrium fluc- 
tuation relations. Left panel: nonlinear conductances for (a) 
Eq. Q, (c) and (e) Eq. (0. Right panel: noise susceptibilities 
for the corresponding relations in the left. Same parameters 
as in Fig. [2] except 74 = 0.1F. 



them as 



kT{^2G aaa -\- G Q _ aQ ), 



(6) 



with G aiaa + G a< aa + 2G aiCc a — 0. In Fig. [3Ja) and (b) 
we explicitly check that these fluctuation relations hold 
despite broken detailed balance. The relations includ- 
ing derivatives of current cross-correlations at different 
conductors, a and /3, read 



7) 



— S a p,a + Saa,fl — kT(Gp, a a + G a ,af) ~ G a , a p), 

with J2 a G 7: a/3 — 0. Eq. ([7]) is verified in Figs. EJc-f). 
It is important to realize here that full access to the fluc- 
tuation relations is only possible in the presence of drag 
current [see Fig.^a) and (c)]. In other words, only when 
detailed balance is broken and a drag current appears are 
all fluctuation relations nontrivially verified. In contrast, 
the absence of drag, i.e. G a ^p = G a ^ = 0, implies 
= 0, for any terminal 7, in which case the fluc- 
tuation relations (JJJ are simply reduced to the relation 

S aa ,f) — 2kTG a ,af3, with Ga.afi — —Ga.ufi- 

Conclusions. — In summary, we have proposed a geom- 
etry of two conductors put in proximity interacting via 
long-range Coulomb forces to test fluctuation relations 
in the non-linear transport regime. This system exhibits 
a drag current as a direct consequence of the absence of 
detailed balance. Our main findings are (i) the general 
expression for the current-voltage characteristic of two 



interacting conductors, and (ii) the verification of the 
fluctuation relations in a nonequilibrium system when 
detailed balance is broken. Our proposal motivates new 
experiments to test the fluctuation relations away from 
equilibrium when detailed balance does not hold. 
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